Viscous Decoupling Transitions For Individually Dragged Particles in Systems with 

Quenched Disorder 
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We show that when an individual particle is dragged through an assembly of other particles in the 
presence of quenched disorder, a viscous decoupling transition occurs between the dragged particle 
and the surrounding particles which is controlled by the quenched disorder. A counterintuitive 
consequence of this transition is that the velocity of the dragged particle can be increased by 
increasing the strength or density of the quenched disorder. The decoupling transition can also 
occur when the external drive on the dragged particle is increased, and is observable as a clear 
signature in the velocity-force response. 

PACS numbers: 82.70.Dd,74.25.Qt 



Recently, several experiments have demonstrated the 
effectiveness of a local probing technique in which an in- 
dividual particle is dragged through an assembly of other 
particles [l|, l^i S Ij] ■ The velocity- force response of the 
driven particle was measured for a single magnetic col- 
loid pulled through a background of nonmagnetic colloids 
in Ref. [1], while in other experiments, optical tweezers 
were used to drive a single colloid through a collection 
of charge stabilized colloids [3|, Ij]. Several numerical 
works have predicted that a rich variety of dynamical 
regimes can occur in these types of systems depending 
on whether the probe particle is driven through a glassy 
media [a, la, [Zl, [sj, |9|, a crystal [lO|, or an assembly of 
rods [ll|. The colloidal dragging experiments are per- 
formed on the microscale; however, it is now becoming 
possible to conduct similar experiments on the meso- and 
nanoscales. For example, recent experiments have shown 
that individual vortices in type-II superconductors can 
be manipulated using a magnetic force microscope tip 
[12, ll3| . This technique could be used to study vortex 
entanglement [l4| and the depinning of vortices from ex- 
tended defects [15[. A simpler initial study would be to 
measure the drag forces on a single driven vortex in the 
presence of other vortices and quenched disorder from the 
sample. Currently, no predictions exist for the response 
in this regime. 

In previous studies of single colloids driven through 
colloidal assemblies, the driven or probe colloid inter- 
acted only with the surrounding colloids. In the super- 
conducting vortex system or in colloidal systems contain- 
ing random quenched disorder, both the probe particle 
and the surrounding particles interact with the underly- 
ing quenched disorder which acts as a pinning potential. 
The addition of pinning might be expected merely to in- 
crease the overall drag on the probe particle; however, 
in this work we show that pinning can induce counter- 
intuitive changes in the probe particle motion since the 
pinning couples to both the probe and the surrounding 
particles. For example, increasing the strength or den- 
sity of pinning sites can induce a decoupling transition 



which sharply reduces the effective drag on the probe 
particle. The same decoupling transition occurs when 
the driving force on the probe particle is increased, and 
appears as a clear jump in the velocity-force curves. In 
the strongly damped regime, the probe particle couples 
strongly to the surrounding particles and induces irre- 
versible plastic deformations. When the pinning strength 
is increased, the surrounding particles are trapped by 
pins and cannot respond to the probe particle, so the 
probe particle decouples from the surrounding particles 
and the effective damping is reduced. Similarly, when 
the drive is sufficiently large, the probe particle is unable 
to induce topological rearrangements in the surrounding 
particles and a decoupling transition occurs. We specifi- 
cally examine a system with screened Coulomb particle- 
particle interactions. Experimentally, this corresponds 
to a two-dimensional colloidal system with quenched dis- 
order, such as in Ref. [l6|. We also study interactions 
appropriate for vortices in type-II superconductors. 

We consider N particles in a two-dimensional sample 
of size L X L with periodic boundary conditions in the 
X and y directions. The overdamped equation of motion 
for particle i is 
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where 77 is the damping constant which we set to 77 = 1. 
The force from particle-particle interactions is F', = 
-J^oE^jVy(ry), where Fq = Z*V(4TCeo), Z* is 
the unit of charge, e is the dielectric constant of the 
medium, and the distance between particles located at 

We use a Yukawa po- 



and Tj is r^ 



= r, - r, 



tential, V{rij) — exp{—Krij)/rij, with screening length 
1/k = 2. The quenched disorder is modeled as Np ran- 
domly placed parabolic pinning traps with pinning force 
Fp = T^k" fpi'^ik/rp)Q{rp - rik)vik, where fp is the pin- 
ning strength, rp = 0.25 is the pinning radius, is the 
Heaviside step function, r^fc = jr^ — r^|, r^ is the location 
of pin A:, and Vik — (r.^ — r^)/rik- The pinning density is 
rip = Np/LF' and the particle density is n = N jl? . We 
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FIG. 1: The velocity-force {v — F) curves for a single driven 
probe particle in a system with Up = 0.85 and fp = 2.0 at 
particle densities n = 0, 0.181, 0.292, 0.375, 0.484, 0.53, and 
0.59 (from top to bottom). Inset: The transition region for 
n = 0.375 (bottom) along with the n = 0.181 curve (top). 



have examined several different system sizes and present 
results for L = 48. For the probe particle, the external 
driving force F^^^ = F^xt'^, and for all remaining parti- 
cles, Fgj.( = 0. We construct a velocity-force curve by 
measuring the time-averaged velocity of the probe parti- 
cle, V = (v ■ x), as a function of Fext- The thermal force 
FJ" arises from random Langevin kicks with the proper- 
ties (Ff) = and (Ff(t)Fj(t')) = 2rikBTS{t - t')S^j. 
Unless otherwise noted, we set T — 0. We also consider 
the case of vortices in type-II superconductors, modeled 
as in Ref. [l7| where F^ 
is a modified Bessel 
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Here 



0o/(27r/ioA^ 



(f) = h/2e is the flux quantum, and A is the London pene- 
tration depth. The probe particle in our system recrosscs 
the same pinned region repeatedly during the course of 
a measurement; to verify that this does not affect the re- 
sults, we tested driving forces applied at low angles to the 
X-axis which cause the probe to cross a different portion 
of the pinned region with each pass through the sample. 
We find the same results for all driving directions. 

We first study a system of Yukawa interacting particles 
and examine velocity force (v — F) relations for the probe 
particle in the absence of thermal fluctuations for varying 
particle densities n. Figure 1 shows the v ~ F curves for 
fixed Up = 0.85 and fp = 2.0 at particle densities ranging 
from n = to n = 0.59. In the limit n = 0, where the 
probe particle interacts only with the pinning, there is an 
initial pinned phase at low drive followed by a transition 
to a moving state at a critical external drive F^^^. = 1.7. 
In the moving state the particle velocity is proportional 
to the damping, v oc ry, representing Ohmic type behav- 
ior. For n < 0.292, the v — F curves are similar to the 
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FIG. 2: Images of a 38 x 12 region of a sample with Up = 
0.85, fp — 2.0, and Fext = 8.0. Large red dot: driven probe 
particle; black dots: unperturbed surrounding particles; open 
circles: surrounding particles perturbed by the passage of the 
probe particle; lines: particle trajectories, (a) The decoupled 
moving regime at n = 0.3 has weak damping, (b) The coupled 
moving regime at n = 0.5 has high damping and significant 
plastic distortions of the surrounding particles. 



n = single particle limit, with only a slight decrease in 
V in the moving state for increasing n. For n > 0.292, in 
addition to the pinned regime at low drive and the Ohmic 
regime at high drive, we find an intermediate non-Ohmic 
regime where the velocity undergoes pronounced fluctu- 
ations and is reduced below the Ohmic value. Between 
the non-Ohmic and Ohmic regimes there is a distinctive 
jump in the velocity at a drive F^^^ indicating a transition 
into a moving state with lower damping and reduced ve- 
locity fluctuations. The inset of Fig. 1 shows a blowup of 
the transition region for n = 0.375; the n — 0.181 curve 
is also presented for comparison. As n increases, the 
velocity of the probe particle in the non-Ohmic regime 
decreases and the transition drive F^^^ increases. 

To illustrate the origin of the different damping in 
the two moving states, in Fig. 2 we show the trajecto- 
ries of the probe particle and the surrounding particles 
in the weakly damped regime at n = 0.3 and in the 
strongly damped regime at n = 0.5. For the weakly 
damped regime in Fig. 2(a), only small distortions of the 
surrounding particles occur as the probe particle moves 
through the sample. In contrast, Fig. 2(b) indicates that 
in the strongly damped regime, the surrounding particles 
are significantly disturbed by the passage of the probe 
particle and undergo irreversible plastic distortion events. 
In order for a plastic rearrangement to occur, a portion of 
the probe particle energy must be transferred to the sur- 
rounding particles, increasing the effective damping on 



the probe particle. The plastic distortions produce large 
fluctuations in the probe particle velocity. In the weakly 
damped regime, the probe particle does not create plastic 
distortions in the surrounding particles. Thus, the transi- 
tion we observe between the weakly and strongly damped 
regimes occurs due to a viscous decoupling transition be- 
tween the probe particle and the surrounding particles. 
The ability of the probe particle to create plastic dis- 
tortions is affected by the external drive, particle den- 
sity, pinning strength, and pinning density. In gen- 
eral, a plastic rearrangement occurs when a surround- 
ing particle is both depinned and displaced by a dis- 
tance ao/2, where qq = ■nT^''^ is the average distance 
between particles. The probe particle spends an aver- 
age time St — a^jF^xt interacting with a surrounding 
particle, and exerts a force F^ = — FoVy(ao/2) on the 
particle. In the absence of pinning, plastic distortions 
would occur when — FoVT^(ao/2)/i^e2^t > 1/2, and thus 
Fl^^ ex -VF(1/(2V^)). The w - i^ curves in Fig. 1 indi- 
cate that a transition between the high damping, plastic 
motion regime and the low damping, decoupled regime 
occurs when the driving force Ff^^t increases. When F^^t 
is large, F^xt > F^^t, the probe particle passes the sur- 
rounding particles so rapidly that it can only induce small 
displacements due to the short interaction time St. These 
displacements are too small to permit plastic rearrange- 
ments to occur, and the probe particle decouples from the 



surrounding particles. For F^^^ < F^. 






plastic 



distortions occur and the probe is in the high damping, 
coupled motion regime. The lower limit of the coupled 
motion regime is determined by the fact that the driving 
force Fext must exceed the critical force F^^t so that the 
probe particle itself remains depinned. 

Figure 1 shows that the transition drive F^^t increases 
with increasing particle density n. This is a result of 
the increasing strength of the interaction Fd between the 
probe particle and the surrounding particles with increas- 
ing n. At low n, the average distance ao between particles 
is large and Fd is low, so the transition to decoupled mo- 
tion Fg^j falls at low values of Ff.xt- As n increases, Fd 
increases due to the decrease in Oo, and the plastic distor- 
tions persist up to higher values of F(.xt: corresponding 
to a higher value of F^^t- In Fig- 3(a) we present a phase 
diagram for F^xt versus n taken from a series oi v — F 
curves. The three phases, pinned (P), decoupled moving 
(DM), and coupled moving (CM), are identified based on 
the values of F^^t ^^^^ Fl^t- The range of driving forces 
over which coupled motion can occur increases with in- 
creasing n. For fixed drive, a velocity drop occurs when 
n is increased above the coupling-decoupling transition, 
as illustrated in the inset of Fig. 3(a) for F^xt = 10. 

To determine the influence of the pinning on the 
coupling-decoupling transition, we vary the pinning force 
fp in a system with fixed n — 0.292 and Up — 0.85, and 
plot the resulting phase diagram in Fig. 3(b). The depin- 
ning force increases linearly with fp, F^^^ ex fp. For very 




FIG. 3: (a) Phase diagram F^xt versus n at fp = 2.0 and 
Up — 0.85 showing the pinned (P), coupled moving (CM), 
and decoupled moving (DM) regimes. Squares: depinning 
threshold F^^f Circles: decoupling threshold F^^^. Inset: 
V versus n at Fext ~ 10.0, showing a drop in v above the 
transition density of n = 0.3. (b) Phase diagram for Fext 
versus fp at Up — 0.85 and n — 0.292. Inset: v versus fp at 
Fext ~ 12.0 shows that v can increase with increasing fp. (c) 
Phase diagram for Fext versus rip for fp — 3.0 and n = 0.292. 
Inset: v versus rip for Fext ~ 12.0. (d) Phase diagram for 
Fext versus T with rip = 0.85, fp = 3.0, and n = 0.292. 
Left inset: ii vs T for Fext = 16.0. Right inset: v vs Fext for 
superconducting vortices with n = 0.44 and rip = 0.46. Lower 
right curve: fp = 0.25; upper right curve: fp = 1.5. 



low fp the system responds as if there were no pinning, 
a coupled motion regime which was previously explored 
in Ref. [5|. At low fp, the probe particle is strongly cou- 
pled to the surrounding particles and generates signifi- 
cant plastic distortions. Since the pinning force tends to 
localize the particles and therefore competes with Fd in 
determining whether the surrounding particles can un- 
dergo plastic motion, we expect plastic distortions to oc- 
cur in the presence of pinning when (Fd — fp)6t > ao/2, 
giving F*xt cc {Fd — fp). This is in agreement with 
Fig. 3(b), where F^^t decreases linearly with fp. The 
coupled motion disappears for fp > 2.0 when the probe 
particle can no longer depin the surrounding particles. 

For fixed drive, the probe particle velocity v initially 
decreases with increasing fp, as shown in the inset of 
Fig. 3(b) for Fext = 12.0. The pinning increases the ef- 
fective viscosity in the coupled moving regime compared 
to the fp = case by increasing the number of irreversible 
events that occur. As fp increases, a surrounding parti- 
cle that has been depinned and displaced by the probe 
particle can more easily be repinned at a new pinning 
site, rather than returning to its previous position in a 
reversible event. When fp is further increased, a sharp in- 
crease in V occurs at the decoupling transition where the 
plastic distortions are lost. We note that v in the decou- 



pled regime is higher than the value of u at /p = 0. These 
results show that increasing the strength of the pinning 
can cause a counterintuitive increase in the velocity (or 
decrease in the damping) of the probe particle. At large 
enough /p, the probe particle itself becomes pinned. 

In Fig. 3(c) we show the phase diagram for F^xt ver- 
sus pinning density Up for a system with fp = 3.0 and 
n — 0.292. There is little change in F^^f with Up, while 
F*^j.f decreases with Up until the coupled motion regime 
disappears for Up > 0.75. When Up increases, the average 
spacing between adjacent pinning sites decreases, and a 
surrounding particle that has been depinned by the probe 
particle can be trapped by a new pinning site before it 
has moved far enough to allow a plastic rearrangement to 
occur. The inset of Fig. 3(c) illustrates the velocity of the 
probe particle versus Up at fixed F^xt = 12.0. For very 
low Up the behavior is similar to the pin-free case. The 
effective viscosity increases sharply at Up « 0.125 when 
the system enters the strongly damped, coupled motion 
regime. For rip > 0.5, the viscosity drops sharply at the 
onset of the decoupled motion regime, and then monoton- 
ically decreases for increasing Up. Increasing the pinning 
density can counterintuitively increase the velocity of the 
probe particle v above the value at np = 0. 

The phase diagram for F^xt vs temperature T is shown 
in Fig. 3(d) for a system with fp = 3.0, Up = 0.85, and 
n = 0.292. At T = the probe particle undergoes de- 
coupled motion for Fext > Fext- ^^ ^ increases, the 
effectiveness of the pinning is reduced, F^^^ decreases, 
and the probe particle recouples to the surrounding par- 
ticles. This produces a velocity drop with increasing T, 
as shown in the left inset of Fig. 3(d) for Fext = 16.0. A 
transition back to the decoupled moving regime occurs 
for Fext > Fext ■ The overall structure of the phase dia- 
gram in Fig. 3(d) is very similar to the dynamical phase 
diagram for vortices in type-II superconductors moving 
over random disorder, where a transition from plastic to 
elastic flow occurs as the driving force on the vortex lat- 
tice is increased 18|. Our results indicate that a similar 



effect can occur even when only a single particle is driven. 

We have performed similar simulations for vortices in 
type-II superconductors. In the right inset of Fig. 3 
we plot the v — F curves from a vortex sample with 
Up = 0.46, Uy = 0.44, and rp = 0.2 in the strong pinning 
regime with fp = 1.5 and in the weak pinning regime with 
fp = 0.25. The v — F curves have the same trend seen in 
Fig. 1 and Fig. 3(b). The fp = 0.25 sample exhibits cou- 
pled motion with high damping, while the fp = 1.5 sam- 
ple is in the low damping, decoupled motion regime. A 
decoupling transition occurs with increasing driving force 
for the weakly pinned sample. This suggests that our re- 
sults should be generic to single driven probe particles 
moving through a background of repulsively interacting 
particles in the presence of quenched disorder. 

In summary, we show that when a single probe par- 
ticle is driven through an assembly of other particles in 



the presence of quenched disorder, a novel viscous de- 
coupling transition can occur between the probe particle 
and the surrounding particles which is controlled by the 
strength and density of the quenched disorder. This tran- 
sition is from a highly damped regime, where the probe 
particle depins the surrounding particles and produces 
irreversible plastic distortions, to a low damping state, 
where the probe particle does not couple to the surround- 
ing particles. Increasing the pinning strength or density 
reduces the coupling between the probe and background 
particles, producing the counterintuitive result that in- 
creasing the strength or density of the quenched disorder 
increases the velocity of the probe particle. The decou- 
pling transition appears as a clear signal in the effective 
damping on the probe particle and produces a distinct 
feature in the velocity force curve. In addition to colloids 
with Yukawa interactions, our results should be general 
to other systems of interacting particles with quenched 
disorder, including vortices in type-II superconductors. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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